Quantum entanglement in random physical states 
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Most states in the Hilbert space are maximally entangled. This fact has proven useful to investigate - among 
other things - the foundations of statistical mechanics. Unfortunately, most states in the Hilbert space of a 
quantum many body system are not physically accessible. We define physical ensembles of states by acting on 
random factorized states by a circuit of length k of random and independent unitaries with local support. This 
simulates an evolution for finite time k generated by a local (time-dependent) Hamiltonian. We apply group 
theoretic methods to study these statistical ensembles. In particular, we study the typicality of entanglement by 
means of the purity of the reduced state. We find that for a time k — 0(1) the typical purity obeys the area law, 
while for a time k ~ O(L) the purity obeys a volume law, with L the linear size of the system. Moreover, we 
show that for large values of k the reduced state becomes very close to the completely mixed state. 



Introduction. — Entanglement is the defining characteristic 
of quantum mechanics: it is a key ingredient in quantum in- 
formation processing! 1 1, quantum many -body theory Q, and 
the description of novel quantum phases of the matter 01 . 
More recently, quantum entanglement has shed new light on 
the foundations of statistical mechanics, and the processes of 
equilibration and thermalization. The idea is simple: even 
though quantum evolution keeps an isolated quantum system 
always in a pure state, the reduced state of a subsystem is 
mixed due to entanglement. One can then speculate whether 
the expectation values of local observables are typically close 
to those of the thermal state [4-<3|. It has been shown that 
such typicality is related to the typicality of entanglement J6). 
For instance, drawing a random state (with the Haar measure) 
from the Hilbert space of a quantum many-body system within 
some energy shell, its reduced state will be typically close to 
the Gibbs state at some temperature T. This is a kinematical 
argument: most states in the Hilbert space look locally ther- 
mal. With a similar reasoning, one can show that starting from 
a non thermal state the infinite time evolution for almost all 
initial states and almost all reasonable Hamiltonians will al- 
most always produce a typical state in the Hilbert space, and 
then the kinematic argument follows. A series of beautiful and 
rigorous results has been found following this approach |6|. 

The problem with this approach is that typical states are 
not physical because they are not accessible in nature. What 
one can access are ground states or low energy states of lo- 
cal Hamiltonians, which are very little entangled because they 
obey an area law [7], and those that can be reached by means 
of evolution with a local Hamiltonian and for a reasonable 
amount of time, say, a time that scales at most polynomially 
in the size of the system. It is clear that the states that can be 
accessed in this way are a set of measure zero in the Hilbert 
space. One needs a doubly exponential time in the system size 
to access all the states of the Hilbert space, which is an unrea- 
sonable time scale for a macroscopic system. For this reason, 
some authors have argued that the Hilbert space is an illusion 
needed for mathematical consistency JSJ. 

On the other hand, physical states do thermalize, as it has 
been shown in several experiments with cold atoms, theoret- 



ical models and numerical simulations ff9l 4TTI . or show typi- 
cality in the expectation value of local observables [ 12]. Does 
this mean that the mechanism for thermalization is not en- 
tanglement? Or the typical entanglement present in physical 
states is still enough to carry over the arguments about ther- 
malization? There are several examples of physical relevance 
in which it is found that when the evolution time scales with 
the size of the system, the state is entangled with a volume law 
|fT3l [T4|. But how typical is this? Can we prove any statement 
about the typicality of such situations? 

In this letter, we want to answer the following question: 
how much are typical physical states entangled? To this 
end we define an ensemble £ ( fc ) f physical states in this 
way: pick a product state of a multi-partite system with lo- 
cal d— dimensional degrees of freedom, then act on this state 
with k independent random unitaries each of them compati- 
ble with some locality structure e.g., supported on edges of a 
graph. These stochastic circuits of local unitaries are, on the 
one hand, mimicking the continuous evolution generated over 
a time A: by a local ( possibly time-dependent) Hamiltonian; 
on the other hand, they are amenable to an elegant analytical 
treatment. In fact, by applying the group theoretic techniques 
of Ref [ 15 1 we are able to compute the ensemble average and 
variance of the purity of a sub region. The result is that, typi- 
cally, for k — 0(1) we obtain an area law, while for k scaling 
as the linear size L of the system the average purity shows a 
volume law. Moreover, we show that fluctuations are small, 
and that there is measure concentration around the average 
value. As a final result, we show that for k — > oo, the subsys- 
tem typically reaches the completely mixed state. 

Statistical ensemble of physical states. — We start by defin- 
ing the ensembles £( k >; henceforth we will refer to the ele- 
ments of these ensembles as the physical states. Let V be a 
set of vertices endowed with a probability measure p : X C 

V f — y p(X) e [0, 1] with J2xP( X ) = L To each of the v^- 
tices x £ V we associate a local d— dimensional Hilbert space 
T~L X ~ C d . The total Hilbert space is thus Hv = ®xevH x , or 
the space of n qudits. A completely factorized state in Hv has 
the form |<f>) = ® x& v \4>x) and let uj = |$) ($| be its density 
matrix. 
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The statistical ensembles of quantum states £^ are con- 
structed in the following way: We first draw a subset X C 
V according to the measure p and then we draw a unitary 
Ux G U(Hx) according to a chosen measure dp(U\X). 
Then we define £(p,dp) = {Ux \^)}x,u x .<s>- This en- 
semble can be generalized to the A;— iterated £^ by con- 
sidering unitaries of the form U — ]X=i Ux t where the 
Xi's (J/xJ's) are drawn according the product probability 

Ui=iP( x i) (IiiLi^C^x,)) : each tick of the clock (la- 
belled by i) a new independent set AQ's and a unitary Uxt 
are picked. In this ensemble, we can compute the statistical 
moments of any hermitian operator. It turns out that by vary- 
ing over the Ux one already picks all the possible factorized 
$'s so that the $ dependence can be in fact dropped. 

Subsystem purity. — We want now inquire about typical en- 
tanglement in the statistical ensembles £^ k \ Let us consider a 
bipartition V = A U B in the system: Hv = Ha ®Hb, 
where Hj = ® x e,jH x with obvious notation. We take a 
state p € £k and consider the reduced density matrix p A — 
Tr b(p)- In order to evaluate the entanglement of p we will 
compute the purity P = Tya(pa)- To compute this trace 
we use the well-known fact that the trace over the square 
of every operator can be computed as the trace of two ten- 
sored copies of that operator times the swap operator. In- 
deed, defining p® 2 = [/® 2 w ® 2 {7t<S>2 an( j considering the or- 
der two shift operator (swap) on T x : H® 2 H> H® 2 , we 
then have P = Tr A (pf 2 f A ) = Tr[p® 2 T A ] where f A = 
® iel T» : (H A )® 2 -> (H A )® 2 by i a > -> \h,h) is 
the order 2 shift operator in the A space Ha — ®i^ A Hi 
and T A _ ■ (H A <8> Hb)® 2 -> (H A ® H B )® 2 is given by 
T A — Ta <8> 1 b ■ We can now consider different concrete en- 
sembles. As a first basic model, let us consider the 

Single edge model. — The system £ edge consists of two sites 
A = {i} and B = {j} connected by an edge e so that the 
Hilbert space is H e = Ha®Hb- The probability distribution 
is the trivial p(e) = 1 and we pick the unitaries U e (H e ) with 
the Haar measure: dp(U e ) = dp Haar- Following Ref. lfTS] 
one can exploit the group theoretic structure of the ensemble 
£ to compute average and statistical moments of operators. 
The average of an operator over a group action is indeed the 
weighted sum of projectors onto the IRreps of the representa- 
tion of that group. The weights are the traces of the operator 
multiplied by the irrep projectors IT, (divided by their dimen- 
sions dj). In our case the IRreps are carried by the totally 
symmetric (H\) and totally antisymmetric (H 2 -) subspaces 
of Hf 2 . The average is given by 

P U = J dUTv [p® 2 T A ] = Tr [u® 2 J dU(U^) m T A U m ] 

Performing the integration over the unitaries we get P U = 
d+ 1 Tr(n+ T)Tr(w 2 n+) = rf+ 1 Tr(n+ T) since uj® 2 is sup- 
ported only in the totally symmetric subspace II + whose di- 
mension is d + = d 2 (d 2 + l)/2 and Tr[a;® 2 II + ] = 1. The 
projector onto the totally symmetric space has the form II + = 
(1 + % Tj)/2 =>> l/2Tr((l (lj) ® 2 + T t T )T t ) = d 3 and 



we finally get P = d 3 /d+ = 2d/(d 2 + 1) = 2N d . In US 

it was defined the average entangling power e p (U) as the 
average entanglement one attains from a factorized bipartite 
state by averaging over the unitaries in the whole space with 

the Haar measure. With this definition, e p (U) := 1 — P U — 
l-2N d = (d- l) 2 /(d 2 + l)[see Eq. (5) in[T5)]. 

Superoperator formulation. — Going back to the general 
case, it is useful to define the superoperator that averages 
over the unitaries U x € U((H A ® H B )® 2 ) K>x{T A ) ■= 
J d/i([/|X)([/^)® 2 T A (C/ x )® 2 - Notice that when dp is the 
Haar measure the TZx are projection superoperators; in the 
rest of the letter we will focus on this case. Then we can eval- 
uate the average purity as 

P=(lo® 2 ,TZ(T a )) (1) 

where 1Z — J^xcv P(X)7Zx is a self-dual (hermitean) super- 
operator. As a far as the purity calculations are concerned this 
superoperator completely characterizes the ensembles £ l - k \ 
Indeed, it is now easy to see that -in view of the statisti- 
cal independence of each iteration- the average purity for 
the k— iterated ensemble £^ is given by the expression ([!]) 
with 1Z replaced by 7Z k . In order to understand the spectral 
properties of 1Z observe that: \\1Z\\ < J2xcv p(X)\\7Zx | < 
J2xcv P( X ) = L Since ^(1) = 1 we then see that \\1Z\\ = 
1 whence the eigenvalues A Q of 7Z are bounded in modulus by 
one and the highest one is Ai = 1, One can then write Pfc = 
T, a ^ k a c a \ where c a := (lu® 2 , (T A ) a ) and (T A ) a denotes the 
projection of Ta onto the eigenvalue X a eigenspace of 7Z. For 
k —> oo this quantity goes to the limite value c\ while the con- 
vergence rate is dictated the second highest eigenvalue A2 of 
K: Pk-Poo ~ c 2 exp(-fclogl/|A 2 |), (k — > 00). 

One of the key steps to obtain the results of this letter is to 
realize that the IZ's superoperators can be regarded as maps 
on the 2^ v > -dimensional space spanned by the Tj^'s (X C V) 
into itself (instead of maps of the tf^l-dimensional C(H^ 2 ) 
into itself). For example, if X :— {a, b} i.e., an edge and A is 
any subset of V, a calculation similar to the above leads to 

1Z X (T A ) = N d (T A \ X + Taux) XnA^$AXnB^(/} 
K X (T A )=T A otherwise (2) 

Moreover, the edge TZx's have low-dimensional invariant 
subspaces of permutations, e.g., in a chain topology the span 
of the T A s associated with connected A's is invariant. This 
remark along with |[T6l allows for drastic simplifications in the 
evaluation of the average purity of £ ( k \ 

To illustrate the usefulness of the algebra |2]), let us compute 
again the purity for the single edge model by means of TZx- 
Here, A is just one site and therefore A\e — % and A U e = 
{A, B}. Moreover, T A (%) = I and 7Z e (T A ) = (I + T A T B ). 
Finally we get P = (uj® 2 ,T a ) = 2N d [Tg). In the case of 
qubits, d = 2 and 2A^ = .8. It is also possible to compute the 
variance by generalizing the group averages to higher power 
of the density operator to obtain AP = .017. A systematic 
treatment of how to compute variances in these statistical en- 
sembles is to be found in IT7l . 
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Random edge model. — At this point, we can implement 
some graph theoretic structure in the set V to make sense of 
the notion of locality for the unitaries U. Consider now the k- 
random edge model £^} ndom defined on a graph T = (V, E). 
We define a flat probability distribution on the edges of the 
graph T: p(X) = l/\E\ if X E E and zero otherwise. Then 
we pick the unitaries on the edges with the Haar measure: 

dfl(U\X) = dflHaar- 

We call DA C E the subset of edges that have non null 
intersection with both A and B. The probability of an edge 
to belong to OA is thus q = \dA\/\E\. We are interested in 
thermodynamical situations where q -C 1. Exploiting Eq. Q 
we see that 



k(t a 



P(X)[T AUX +T A \ X ] 



E pW t * + E 

XeE\dA XedA 

(3) 

where X g E is an edge of the graph. For k = 1 is then easy 
to find the average purity: P = Y^xeE\dAP( X )( u ® 2 > T a) + 
ExedA N d (w® 2 ,T Aux + T A \x) = 1 - (1 - 2N d )q. From 
Eq.Q we see that every application of 1Z X transforms the 
subset A into a superposition of AUX and A\X so that at any 
successive iteration the boundary of the new subset changes 
and its boundary length may change. To understand the struc- 
ture of 1Z k (Ta) it is instructive to consider explicitly the k = 2 
case. Iterating Q and computing the purity using ([T]i [116] we 
find P ~ (1 - q) 2 + (1 - q)2qN d + (2qN d ) 2 = (1 - q(l - 
2N d )) 2 where we have used \d(A U e)| ~ \d(A\e)\ ~ \dA\. 
This is true e.g. if every vertex has degree o(|cM|). The cal- 
culation easily extends to the case k = o(\dA\) and one finds 
P~k ~ (1 - q(l - 2N d )) k . Therefore the average 2-Renyi 
entropy of the ensemble is 82(h) 
-k\og(l-q(l-2N d )) 



logP fc > -log'Pfe ~ 

(1 — 2N d )qk. In terms of the aver- 

V 



age entangling power one gets 82(h) > k\dA\/\E\ e p (U) 
In words: the average 2— Renyi entropy for the random edge 
model of the fc— th iteration is lower bounded by k times the 
fraction of vertices in the boundary of the region A times the 
average entangling power of an edge unitary. This shows that 
there is a linear increase of entropy in time. Moreover, one can 
compute variances of P and show that AP ~ \dA\ ifTTl so 
that v AP/P ~ l/y/\dA\. This in turn implies measure con- 
centration (typicality) in the thermodynamic limit \dA\ — > oo. 
In the extreme case in which the graph is complete, then q = 1 
and one would expect that for k large the completely mixed 
state is reached [18|. 

This random edge model is of course ultra-local and then 
somewhat unphysical. The associated ensemble, roughly 
speaking, corresponds to the set of states obtainable from a 
product state by the action of a time dependent Hamiltonians 
supported on a single edge. We now move to consider a more 
physically satisfactory scenario. 

The linear chain. — We now consider a model in which all 
the unitaries act on all the edges of the graph T. The prob- 
ability distribution is thus p(X) = 1 for X = V and zero 
otherwise. For the sake of simplicity in the following, we 
will consider the case of the graph V being a bipartite chain 




FIG. 1. A bipartite (A,B) spin chain of length L = La + Lb 
with nearest-neighbor qubits interacting via 2— qubit gates (ellypses). 
The edge e is the one that straddles the two partitions. The gates 
are numbered by the subscript xi where x — A,B denotes the two 
halves of the chain and i the distance from the boundary of the two 
partitions. 



of length L = La + Lb ■ Extensions to higher dimensional 
geometries will be presented in ifTTlO . We will label by U e 
the unitary acting on the edge straddling the (A, B) biparti- 
tion, while we will use the labels ai , bi for the unitaries that 
act in the bulk of A, B respectively (see FigJTJ. We label the 
sites of the chain as La, ..,1aAb, ■■■Lb- Since the unitary 
is a product of all the edge unitaries, we need to specify in 
which order they act. In the following, the unitary U a will al- 
ways denote the product over all the edges in E with the order 
given by the permutation a, so U a is the ordered product over 
local two-qudit unitaries. This corresponds to the (time or- 
dered) infinitesimal evolution with a local Hamiltonian, where 



a gives the time ordering: U a = U c 



cr(ei) • ■ • ^ CT (e| E |) 



At this 



point we construct the set fS(T) = {U a \&)}u with mea- 
sure dfi(U) = S(U — Ua) Ilees d^Haar(U e )- This ensemble 
approximates all the states that can be evolved from a factor- 
ized state with a local Hamiltonian acting for an infinitesimal 
amount of time. By fc-iteration, we obtain the time evolution 

for a finite time k : £ ( ^ am (T) = {nti ^ l*> W- Here we 
consider all the possible ordered sequences of unitaries by tak- 
ing, a each time step, a permutation a of the edges uniformly 
at random. 

This ensemble approximates all the states that can be 
reached in time k by the evolutions originated by all the pos- 
sible time-dependent Hamiltonians on a graph. The ensem- 
ble £ thus only depends on the number of iterations (the 
"time") k and the graph V. The loss of purity due to the ac- 
tion of the unitaries depends on their order. Thus, in order 
to find an upper bound to the average purity, we consider the 
ordering that gives the minimum loss of purity. We can see 
that, for k = 1, acting in A and B after having acted on the 
edge e does not change the purity, so a sequence of the type 
UaUbU e results in a minimal loss of purity. Moreover, the 
order of the unitaries inside A and B also counts. From the 
iteration of the algebra Eq.|2} we can see that if we pick the 
ordering in which we first act near the boundary and proceed 
towards the outer parts of the chain: UaUbU c where Ua = 
Ua L U aL ...U a2 U ai and U B = U b U b ....U b2 U bl 
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we will get the lowest possible powers of Nd and correspond- 
ingly the least decrease of purity. As k increases, the differ- 
ence between different orderings is attenuated and for very 
large values of k it can also be neglected. We will anyway 
always consider the worst case scenario of ordering a which 
corresponds to the minimal decrease of purity. 

The action of the superoperator 1Z in the present chain 
case is more complicated because now X is not just the sup- 
port of one unitary, but it contains the ordered product of 
all the edges. In particular, notice that now 7Z is not her- 
mitean. Using Eq.Q multiple times we find: for k = 1, 
H(T A ) = 7V d (TA_i + T A+1 ). Where we used the notation 
A + r = AU{l B ,...,r B }and A-r = A\{1 A , .., r A }. At 
the second iteration k = 2, we get 

K 2 (T A ) = N 2 T A _ 2 + 2NlT A ^ + 2N*T A+1 + N 2 d T A+2 . 

We can see that nodes at distance 2 from the boundary enter 
the expression. Each T in the expression for lZ k (T A ) gives 
a 1 when we take the scalar product with w® 2 . So we find 
P fc=2 = 2iV| + 4iVj. It is important to understand how the 
interactions propagate with k. A somewhat lengthy calcula- 
tion shows that as k increases, nodes at distance k from the 
edge participate to the averaging procedure and for every node 
that participates we pick a power for the base Ng, The result 

of the calculation is that P k = Y.ZlV H^IT ^ fc+m) . 
This is valid for k < L A . Summing the series for P k for large 
values of k one finds P k ~ 2 [N d /(1 ~ N d )] k . Recall that 
this equation, in view of the choice of the a corresponding 
to the U u with the least entangling power, is an upper bound 

(k) 

for the purity in £ chain - The exponential decay of the purity 
in k is due to the fact that all the qudits at distance k from 
the edge are getting mixed. The average 2— Renyi entropy is 
52 >= — log Pk, therefore we have the bound 

S 2 > k log - log 2 ~ fc log d ~ log 2 (4) 

(Last approximation holds for large d). Eq. (|4]) for k = 
0(L A ) implies a) a volume law for the entanglement scaling 
b) typicality: a nearly minimal value of the average of purity 
(in view of the Markov inequality), forces also the fluctua- 
tions around this average to be small. For k > \L A \ one has 
no longer a linear increase of entanglement with time but ob- 
serves a saturation. This is the type of behavior that has been 
found in examples of entanglement dynamics after a quench 
using CFT techniques |[T3l[T9l [see also the recent papers |20l . 
closer to the spirit of the present one]. 

To study the limit of average purity for k —> oo we first 
notice that the chain superoperator is a (er-ordered) product of 
(non-commuting) projections TZ cha i n = *R><t(\e\) ' ' ' ^o-(i) ■ 
This implies \\TZ chain \\ < YieeE \\"^e\\ < 1; again this means 
that all the eigenvalues of lZ c hai n are smaller in modulus than 
1 and therefore asymptotically just fixed points e.g. ,1,7V, 
contribution to T A survives. If now one assumes that the sym- 
metric combination 1 +Ty is the only relevant fixed point one 
finds Pfe^oo = (d 2L - hA +d L+LA )/d L {d L + 1). We have 



checked this result by numerical simulations [17] for the least 
and most entangling er's but we conjecture it to hold true for 
all orderings and besides the one-dimensional chain scenario. 
For large \V\ = L one has P k ^oo — d~ LA + d hA ~ L that in 
turn for L A < L/2 shows that the asymptotic purity differs 
from that of the totally mixed state t A jd hA for terms of order 
dr LA . Finally, if L A ^> 1, this implies that the vast majority 
of the states in E^ k \ once reduced to A, are close in Li-norm 
to the maximally mixed state. 

Conclusions.-. In this letter we addressed the problem of 
typical entanglement of physical quantum states. This issue 
is crucial to assess the relevance of arguments based on en- 
tanglement typicality for generic random states in the foun- 
dations of statistical mechanics and equilibration dynamics. 
To this aim we have defined statistical ensembles of physi- 
cal states by considering product states on a multi -partite 
systems and evolving them with k independent stochastic lo- 
cal gates. Ensemble averages can be computed by introducing 
suitable superoperators and using group-theoretic tools as in 
lfT31l . For one-dimensional geometries (chains of length L) we 
proved that a) For a finite k = 0(1) time, the entanglement 
measured by the decrease in purity shows an area law., b) As 
the time increases, correlations propagate and for k = O(L) 
the size of the system states display a volume law. We ex- 
pect these results to hold true for higher dimensional regular 
graphs e,.g., 2D lattices lfl7l . We would like also to stress that 
even if in this Letter we have used purity to quantify the en- 
tanglement, this method extends in a straightforward way to 
general a-Renyi entropy by natural modifications of superop- 
erator 1Z and permutation T A in Eq. <|TJ [ 1 7 1 . As we have seen, 
statistical fluctuations are small: the entanglement behaviour 
described above is typical and typicality of entanglement is 
at the root of the local thermalization results. At this point, 
we may speculate as whether our results implies a local ther- 
malization result for physical states. This would be important 
both for the thermalization problem and for quantum informa- 
tion processing. 
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